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A  high  strength  magnesium  alloy,  AMX602  (Mg-6%Al-0.5%Mn-2%Ca),  was  manufactured  by  the  spinning 
water  atomization  process  (SWAP)  and  extruded  into  bar  and  plate  geometries.  Microstructural  analysis 
using  electron  backscatter  diffraction  revealed  that  the  processing  produces  an  alloy  with  grains  between 
0.5  and  5  pm,  with  comparatively  weak  texture  for  Mg.  The  plate  and  bar  had  different  textures-the 
former  had  a  conventional  hexagonal-close-packed  (HCP)  rolling  texture  and  the  latter  a  HCP  extrusion 
texture.  Quasi-static  and  dynamic  compression  experiments  were  performed  to  probe  the  material’s 
mechanical  behavior  in  the  three  processing  directions.  The  experiments  on  each  geometry  revealed 
different  anisotropic  properties  induced  by  a  change  in  the  active  deformation  mechanisms.  The  ani¬ 
sotropy  was  more  pronounced  at  dynamic  strain  rates  than  quasi-static.  A  reduced-order  crystal  plas¬ 
ticity  model  that  demarcates  twinning,  basal  slip,  and  non-basal  slip  mechanisms  was  fit  to  the  ex¬ 
perimental  data  from  the  plate  and  bar.  The  model  was  consistent  with  experimental  data  and  revealed 
that  in  the  plate  twinning  dominated  yielding  in  the  extrusion  and  transverse  directions,  but  slip 
dominated  the  normal  direction.  Yielding  in  the  bar  was  dominated  by  twinning  in  the  extrusion  di¬ 
rection,  but  both  slip  and  twinning  were  significant  in  the  other  two  directions.  The  model  showed  the 
different  anisotropic  responses  were  due  to  the  different  textures  produced  during  the  processing  of  each 
geometry.  Lastly,  our  data  provided  the  basis  for  considering  twinning  to  be  rate  insensitive  in  the  model, 
which  we  confirm  to  be  valid  to  at  least  5000  s-1. 

Published  by  Elsevier  B.V. 


1.  Introduction 

Magnesium  is  the  lightest  structural  metal  with  a  specific 
strength  similar  to  or  exceeding  that  of  aluminum  and  titanium 
[1  ].  However,  commercial  magnesium  alloys  have  limited  ductility, 
anisotropic  mechanical  behavior,  and  relatively  poor  corrosion 
resistance.  The  strength  and  ductility  can  be  optimized  through 
alloy  development,  texture  modification,  grain  refinement,  or 
combinations  thereof  [2,3],  For  example,  recent  Mg  alloy  devel¬ 
opments  have  consisted  of  enhancing  precipitate  hardening  via 
microalloying  by  adding  trace  amounts  of  certain  elements  to 
create  finely  dispersed  precipitates  [4],  In  addition,  alloys  with  a 
more  random  texture  have  shown  increased  ductility  than  typi¬ 
cally  textured  wrought  alloys  because  the  random  texture  in¬ 
creases  the  strain  hardening  rate,  which  promotes  stable  plastic 
deformation  [5-7],  Furthermore,  reducing  the  grain  size  increases 
the  strength  of  a  metal  through  the  Hall-Petch  relationship. 
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Twinning  also  follows  this  relationship  [8],  but  the  critical  stress  to 
activate  twinning  increases  faster  than  for  slip  as  the  grain  size 
decreases,  resulting  in  twinning  being  suppressed  below  a  critical 
value  [8,9  .  Typical  processing  methods,  e.g.  rolling  and  extrusion, 
limit  the  ability  to  reduce  the  grain  size  due  to  Mg's  limited  duc¬ 
tility,  and  these  methods  produce  sharp  textures  associated  with 
increasing  mechanical  anisotropy,  tension/compression  asym¬ 
metry,  and  limited  ductility. 

Severe  plastic  deformation  processing  methods  are  effective  at 
refining  grain  size,  but  they  also  produce  sharp  textures  [10-13], 
Equal  channel  angular  pressing  (ECAP)  is  able  to  reduce  the  grain 
size  and  tailor  the  texture  for  property  improvement.  Several  au¬ 
thors  [14-16]  have  performed  ECAP  on  AZ31B  which  showed  in¬ 
creased  strength  and  reduced  anisotropy  while  preserving  ducti¬ 
lity  using  specific  processing  routes  and  temperatures.  The  re¬ 
sulting  basal  texture  is  generally  stronger  than  before,  and  is  ro¬ 
tated  with  respect  to  a  principal  processing  direction.  Therefore, 
when  determining  the  subsequent  strength  and  ductility,  the 
loading  directions  were  at  an  angle  with  respect  to  the  basal 
texture.  When  conventional  rolled  AZ31B  is  loaded  along  the 
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processing  directions,  the  loading  direction  coincides  with  the 
principal  texture  components,  maximizing  the  measured  the  ani¬ 
sotropic  response.  The  reduced  anisotropic  behavior  and  negligible 
change  in  ductility  of  ECAPed  AZ31B  is  partially  a  consequence  of 
the  chosen  loading  directions,  and  not  an  intrinsic  material 
property.  When  samples  are  cut  along  principle  material  direc¬ 
tions,  the  measured  yield  strength  differences  in  different  direc¬ 
tions  can  still  vary  by  a  factor  of  2  or  more,  and  limited  ductility  is 
observed  for  particular  loading  conditions  14,15],  In  applications 
with  complex  loading  conditions,  where  the  full  anisotropic  ma¬ 
terial  response  is  probed,  it  is  unclear  whether  this  ECAPed  AZ31 B 
will  exhibit  its  measured  strength  and  ductility  improvement. 

Rapid  solidification  is  another  technique  for  concurrent  grain 
refinement  and  texture  modification.  Several  authors  17—20]  have 
used  rapid  solidification  processes  to  prepare  metal  powders  such 
as  Fe,  Mg,  Al  and  Cu  and  their  alloys.  Additionally,  the  Spinning 
Water  Atomization  Process  (SWAP)  is  capable  of  producing  metal 
powders  with  a  fine-grained  microstructure  [21-24],  Recently, 
SWAP  was  applied  to  Mg  alloys  and  it  has  shown  promise  in 
producing  relatively  large  volumes  with  a  fine-grained  micro¬ 
structure  and  a  comparatively  weak  texture  [25-28],  SWAP  con¬ 
sists  of  rapidly  cooling  molten  material  with  water,  which  pro¬ 
duces  a  coarse  powder.  The  powder  is  consolidated  and  extruded 
to  its  bulk  shape.  Two  candidate  high  strength  Mg  alloys,  AMX602 
(Mg-6%Al-0.5%Mn-2%Ca)  and  ZAXE1711  (Mg-l%Zn-7%Al-l%Ca-l% 
La),  have  exhibited  improved  ballistic  performance  compared  with 
conventionally  processed  AZ31B  [29], 

An  initial  batch  of  AMX602  specimens  was  produced  in  ex¬ 
truded  bars  with  a  width-to-thickness  ratio  of  1.5:1.  These  speci¬ 
mens  had  reduced  peak  texture  (compared  to  conventionally  rol¬ 
led  Mg  alloys),  measured  yield  strengths  of  approximately  250- 
350  MPa,  and  strain-to-failures  of  about  10-20%,  depending  on  the 
loading  direction,  strain  rate,  and  sense  (tension  or  compression) 
[26,27,30].  Because  commercial  applications  need  plate-like  spe¬ 
cimens  with  larger  width-to-thickness  ratios,  a  subsequent  batch 
was  produced  with  a  ratio  of  6:1. 

In  this  study  AMX602  was  prepared  by  SWAP  and  the  coarse 
powder  was  consolidated  and  extruded  into  a  bar  or  plate  with  a 
width-to-thickness  ratio  of  1.5:1  and  6:1,  respectively.  We  ex¬ 
amined  the  microstructure  of  the  plate  using  electron  backscatter 
diffraction  (EBSD).  We  performed  experiments  in  both  geometries 
to  determine  the  material’s  mechanical  behavior  in  compression  at 
quasi-static  and  dynamic  strain  rates  in  each  of  the  three  proces¬ 
sing  directions-extrusion  (ED),  transverse  (TD)  and  normal  (ND) 
directions.  Additionally,  a  recently  developed  HCP  crystal  plasticity 
model  [31  ]  was  used  to  analyze  the  material's  response  subjected 
to  the  aforementioned  loading  conditions.  The  model  was  used  to 


show  that  differences  in  the  mechanical  response  and  anisotropy 
of  the  bar  and  plate  can  be  attributed  to  the  variations  in  the 
texture  between  the  geometries.  We  also  show  experimentally 
that  twinning  is  strain  rate  insensitive  over  the  wide  range  of 
strain  rates  tested  here,  up  to  5000  s  which  is  an  assumption  of 
the  model. 


2.  Experimental  procedures 

2.3.  Material  processing 

AMX602  coarse  powder  was  prepared  via  SWAP,  as  shown 
schematically  in  Fig.  1.  The  entire  process  is  performed  in  an  inert 
atmosphere.  The  molten  AMX602  is  directly  streamed  into  the 
chamber  where  water  is  spun  around  the  interior  surface.  Im¬ 
pinging  water  droplets  rapidly  cool  the  molten  stream,  which 
locks  in  a  fine-grained  microstructure.  The  powder  is  consolidated 
into  a  green  compact  and  extruded  into  the  final  shape.  The  bar 
geometry  investigated  here  is  the  same  as  the  initial  batch  pro¬ 
duced  previously  [26,27,30],  with  a  compaction  and  extrusion 
temperature  of  250  °C.  For  the  plate  geometry,  the  powder  was 
consolidated  and  extruded  at  370  °C  into  a  cross  section  of 
152  mm  x  25.4  mm  (6:1  width-to-thickness  ratio).  There  is  some 
oxidation  of  the  molten  Mg  when  it  interacts  with  the  water,  ac¬ 
cording  to  Kondoh  et  al.  [25],  however  the  concentration  was  not 
quantified.  We  scanned  an  approximately  100  x  100  pm2  area  of 
the  plate  geometry  using  energy  dispersive  x-ray  spectroscopy 
(EDS)  and  measured  an  oxygen  concentration  of  0.8  wt%.  The 
oxygen  EDS  map  (not  shown)  indicated  a  generally  uniform  dis¬ 
tribution  over  the  area  scanned,  which  is  consistent  with  [25  .  In 
this  paper,  AMX602  refers  to  the  final  extruded  material.  The  re¬ 
sponse  of  the  compact  and  extruded  powder  were  not  considered 
in  this  work,  but  was  done  by  previous  authors  [25,27,29], 

2.2.  Electron  backscatter  diffraction  analysis 

A  FEI  Nova  NanoSEM  600  SEM  (FEI  Co.  Hillsboro,  OR,  USA)  was 
used  to  characterize  the  microstructure  of  AMX602.  Electron 
backscatter  diffraction  (EBSD)  characterization  was  used  to  quan¬ 
tify  crystallographic  orientations  in  the  samples  at  a  20  kV  accel¬ 
erating  voltage  and  at  70°  tilt  with  an  EDAX  EBSD  system  (EDAX 
Inc.  Mahwah,  NJ)  in  the  SEM.  EBSD  patterns  were  collected  with  a 
step  size  of  0.2  pm  on  the  ED-TD  plane  in  the  middle  of  the  plate, 
and  approximately  in  the  middle  of  the  plate  thickness  (i.e.  middle 
of  the  normal  direction  of  the  plate).  The  collected  data  was 
minimally  "cleaned"  by  using  the  TSL  OIM  Analysis  7  "cleanup” 
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Fig.  1.  Schematic  of  SWAP  equipment  to  produce  AMX602  powder  (left)  and  morphology  of  the  powder  [25]. 
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program  to  correct  inaccurately  indexed  points  based  on  neigh¬ 
boring  point’s  orientation  correlation.  A  minimum  of  150,000 
grains  are  used  to  quantify  the  texture. 

2.3.  Quasi-static  and  dynamic  compression  experiments 

Samples  cut  from  the  bar  and  plate  were  loaded  to  failure  in 
compression  at  quasi-static  and  dynamic  loading  rates  along  the 
extrusion  (ED),  transverse  (TD)  and  normal  directions  (ND).  Quasi¬ 
static  compression  experiments  were  performed  using  a  screw- 
driven  Instron  machine  at  a  constant  displacement  rate  corre¬ 
sponding  to  a  nominal  strain  rate  of  10~4s~\  Two  strain  gages 
mounted  on  orthogonal  surfaces  measured  strain  on  the  speci¬ 
mens,  which  were  8.5  mm  x  8.5  mm  x  12.7  mm.  Dynamic  com¬ 
pression  experiments  were  performed  on  11.4  and  6.35  mm  cubic 
samples  for  the  plate  and  bar  geometries,  respectively,  using  a 
split-Hopkinson  pressure  bar  (SHPB)  with  19  mm  diameter 
maraging  steel  bars.  The  stress-strain  history  of  the  sample  was 
calculated  from  the  incident,  reflected  and  transmitted  strain  gage 
signals  mounted  on  the  incident  and  transmitted  bars  [32],  The 
nominal  strain  rate  was  approximately  1000  s_1.  All  samples  were 
cut  from  the  center  of  the  plate  and  bar  to  eliminate  any  surface 
effects. 


3.  Modeling  framework 

We  used  a  combined  crystal-  and  macro-plasticity  framework 
to  model  the  quasi-static  and  dynamic  response  of  AMX602  spe¬ 
cimens.  Elsewhere  Becker  and  Lloyd  [31]  describe  details  of  the 
model,  so  herein  we  restrict  ourselves  to  essential  model  features 
and  assumptions. 

The  model  uses  a  pseudo-slip  description  for  extension  twin¬ 
ning  that  neglects  twin  reorientation  effects  on  the  mechanical 
response,  and  assumes  slip-based  mechanisms  are  basal  <a)  slip, 
prismatic  <a)  slip,  and  pyramidal  <c+a)  slip.  Following  Graff  et  al. 
[33]  the  inelastic  velocity  gradient  is  partitioned  according  to  the 
individual  mechanisms,  i.e., 

6  3 

LP  =  2  Cn(s“win  ®  m“win)  +  X  /b“as(Stas  ®  “bas) 
a=  1  a=  1 

3  6 

+  X  Jpris(Spris  ®  mpris)  +  X  /pyrm(Spyrm  ®  mpyrm)' 

0=1  o=l 

where,  for  example,  s“win  and  m"win  denote  the  deformation  plane 
direction  and  normal  of  twin  system  a,  respectively.  The  reduced- 
order  model  is  based  on  simplifications  of  the  above  relations  for 
computational  efficiency  and  stability,  so  that  the  framework  can 
be  used  to  simulate  large  dynamic  problems  of  interest.  The  model 
approximates  the  combined  effect  of  basal  (a)  slip  by  assuming 
basal  slip  occurs  on  one  aggregate  system  in  the  direction  of 
maximum  shear  stress  on  the  basal  plane,  denoted  sbas,  such  that 

3 

2  ^bas(Sbas  ®  mbas)  /bas(  Sbas  ®  mbas)’ 

0  =  1 

where  the  basal  slip  normal  is  identical  for  all  3  systems.  Because 
prismatic  and  pyramidal  have  similar  strengths,  and  their  com¬ 
bined  effect  forms  an  approximately  closed  yield  surface  in  biaxial 
stress  space  [34],  they  are  approximated  using  as  a  single  isotropic, 
rate-dependent  plastic  deformation  mechanism,  i.e., 

3  6  „  , 

X  ipris(Spris  ®  mpris)  +  ^  /pyrm(Spyrm  ®  mpyrm)  ^slip^T;  ■ 
o=l  o=l 

where  eslip  denotes  the  equivalent  plastic  slip  rate,  a'  denotes  the 


deviatoric  portion  of  the  Cauchy  stress,  and  is  the  von  Mises 
equivalent  stress.  Following  these  approximations,  the  plastic  ve¬ 
locity  gradient  for  the  reduced-order  model  is  expressed  as 

LP  =  2  CnKwin  ®  mtwin)  +  /bas(  *bas  ®  mbas)  +  eslipy^-. 
o=l 

Despite  the  seemingly  restrictive  set  of  kinematic  assumptions 
employed  in  the  reduced-order  model,  it  has  been  shown  to  re¬ 
liably  capture  the  mechanical  response  of  single-crystal  Mg 

[31,35], 

Differing  constitutive  relations  and  material  parameters  are 
used  to  describe  the  disparate  responses  of  the  three  deformation 
mechanisms.  Basal  slip  is  assumed  to  be  rate-independent,  and  its 
flow  stress  evolves  according  to  a  Voce-type  hardening  law  [36], 
where 


Tbas  ~  A),bas  +  rbas 


1  -  exp 


^basfbas  j 

Tb"s  j 


Extension  twinning  is  also  assumed  to  be  rate-independent, 
and  the  flow  stress  on  each  twin  system  evolves  as 

Awin  —  A), twin  ^  twin /twin’ 


where  /twin  =  iLiCin'  and  Awin=°°  when  yt“>0.1289.  This  model 
assumes  twin  resistance  increases  linearly  and  proportionally  to 
the  total  amount  of  twinning  that  has  occurred  on  all  twin  sys¬ 
tems,  which  is  equivalent  to  assuming  isotropic  cross  hardening. 
The  resistance  of  prismatic  and  pyramidal  slip  is  assumed  to  be 
rate-dependent,  and  influenced  by  the  extent  of  basal  slip  and 
twinning,  such  that 


°slip  -  |ffO,slip  +  °sBp 

slip 

£o  J 


£slip(^bas/bas  Wtwin)j( ' 


1-exp 


^slip(eslip  +  <Jba/bas) 
^00 
°slip 


In  the  above  equation  gs|ip  is  a  linear  hardening  term  that 
causes  basal  slip  and  twinning  to  cross-harden  the  combined 
prismatic  and  pyramidal  slip  systems,  where  qbas  and  qtwin  control 
the  rate  of  cross  hardening.  Note  that  basal  slip  is  assumed  to  also 
affect  the  rate  that  slip  approaches  its  saturation  stress  in  the 
Voce-type  relation,  i.e.,  <r0  slip  +  crsbp.  Material  parameters  for  the 
constitutive  model  are  given  in  Table  1. 


Table  1 

Material  parameters  for  re¬ 
duced-order  crystal  model  of 
AMX602. 


Parameter 

Value 

Units 

T0,bas 

80 

MPa 

Tbas 

10 

MPa 

^bas 

5*103 

MPa 

T0,twin 

80 

MPa 

^twin 

50 

MPa 

^O.slip 

160 

MPa 

°slip 

250 

MPa 

^slip 

to4 

MPa 

£slip 

125 

MPa 

%as 

2 

- 

^twin 

15 

- 

1 

s-> 

m 

0.01 

- 
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4.  Results  and  discussion 

4.1.  As-received  AMX602  microstructures 

Fig.  2  shows  the  microstructure  of  the  AMX602  plate  specimen 
revealed  using  EBSD.  The  scan  shows  nearly  equiaxed  grains  with 
diameters  ranging  from  0.5  to  5  |im.  Elongated  grains  larger  than 
10  pm  in  size  were  sparsely  distributed  throughout  in  other  scans 
(not  shown).  During  SWAP,  the  rapid  cooling  of  the  molten  ma¬ 
terial  likely  created  the  elongated  grains,  which  solidified  slower 
than  the  surrounding  region  allowing  those  particular  grains  to 
grow  larger.  Additionally,  no  obvious  porosity  is  observed.  The  bar 
geometry  microstructure  is  characterized  elsewhere  [26],  It  has  a 
slightly  finer  grain  size  when  extruded  at  300  °C  and  a  much  larger 
grain  size  when  extruded  at  400  °C. 

We  sampled  random  orientations  from  the  EBSD  maps  ob¬ 
tained  in  this  work  for  the  plate,  and  in  prior  work  for  the  bar 
geometry  [26],  and  these  orientations  are  projected  onto  the  ex¬ 
perimental  {  0002}  and  { 1010}  pole  figures  (Fig.  3).  Note  the  gray¬ 
scale  shading  is  the  same  for  all  the  pole  figures  but  the  scale 
values  are  different  for  each.  In  the  plate  (Fig.  3(a)),  the  basal 
texture  is  orientated  primarily  in  the  ND,  which  is  consistent  with 
a  rolled  HCP  material,  albeit  with  a  much  weaker  peak  intensity. 
The  prismatic  planes  are  preferentially  oriented  toward  the  ED-TD 
plane  and  are  axisymmetric  about  the  ND.  The  basal  texture  of  the 
bar  (Fig.  3(b))  is  essentially  axisymmetric  with  respect  to  the  ED 
and  is  distributed  evenly  along  the  TD-ND  plane.  This  is  typical 
extrusion  texture  for  HCP  metals.  The  prismatic  planes  are  ran¬ 
domly  orientated  with  a  slight  concentration  near  the  ED.  The 
plate  has  slightly  greater  basal  texture  than  the  bar  (4.2  versus  2.9, 
respectively).  The  difference  in  the  textures  between  plate  and  bar 
specimens  is  a  consequence  of  the  extrusion  parameters.  The  in¬ 
itial  compact  for  the  plate  was  280  mm  in  diameter  and  was 
subsequently  extruded  into  a  plate  with  dimensions  of 
152  mm  x  25.4  mm,  which  is  a  reduction  during  extrusion  pro¬ 
cessing  in  the  TD  of  1.8:1  and  11:1  in  the  ND.  The  plate  is  domi¬ 
nated  by  the  reduction  in  thickness  so  the  texture  resembles  that 
of  a  rolled  plate.  In  contrast,  the  compact  for  the  bar  was  reduced 


Fig.  2.  EBSD  map  of  the  AMX602  plate  specimen. 


by  a  similar  magnitude  in  both  the  thickness  and  width  directions, 
which  is  similar  to  rod  extrusion  texture.  This  implies  there  will  be 
a  greater  volume  of  material  with  a  rolled  texture  as  the  width-to- 
thickness  ratio  is  increased  of  the  extruded  geometry.  Under¬ 
standing  how  the  processing  history  influences  the  overall  texture, 
and  therefore  the  material  response,  is  essential  for  designing 
larger  plates  used  in  engineering-scale  applications. 

4.2.  Mechanical  behavior  of  AMX602 

We  performed  compression  experiments  in  the  plate  and  bar  at 
10~4s_1  and  1000  s_1  in  the  ED,  TD  and  ND  to  determine  the 
influence  of  loading  rate  on  the  anisotropy  of  the  mechanical  re¬ 
sponse.  All  experiments  resulted  in  the  macroscopic  failure  of  the 
sample,  which  is  reflected  in  the  stress-strain  curves  as  a  quick  or 
sudden  drop  in  the  flow  stress  at  large  strains.  Multiple  experi¬ 
ments  were  performed  for  each  geometry,  loading  direction,  and 
strain  rate  and  showed  consistent  results,  thus,  a  single  re¬ 
presentative  test  is  shown  for  clarity. 

Fig.  4  shows  the  true  stress-strain  response  for  the  three 
loading  directions  of  the  plate  at  quasi-static  and  dynamic  strain 
rates.  This  figure  also  shows  the  model  simulations  but  these  re¬ 
sults  will  be  discussed  in  the  next  section.  Compression  along  the 
ED  results  in  a  characteristic  s-shaped  hardening  curve,  which  is 
indicative  of  { 1012}  extension  twinning  [371.  For  both  strain  rates, 
after  yielding  the  material  exhibits  a  nearly  perfectly  plastic  stress 
plateau  until  strains  of  approximately  4%.  Additionally,  despite 
seven  orders  of  magnitude  difference  in  the  strain  rate,  there  is  no 
difference  in  the  yield  strength  (~280  MPa)  and  little  difference  in 
the  post-yield  stress-strain  response.  The  dynamic  strain  rate  has  a 
strain  to  failure  that  is  0.04  more  than  for  the  quasi-static  strain 
rate.  Compression  along  the  TD,  Fig.  4(b),  exhibits  a  plateau  region 
that  is  less  pronounced  with  slight  hardening,  and  little  appreci¬ 
able  difference  in  the  yield  strength  (~  240  MPa)  between  high 
and  low  strain  rates.  The  subsequent  work  hardening  rate  is 
slightly  reduced  at  the  quasi-static  strain  rate  and  further  reduced 
as  compared  to  the  ED.  Therefore,  although  the  strains  to  failure 
are  similar  for  compression  along  the  ED  and  TD,  the  peak  stress  at 
failure  is  lower  along  the  TD  than  the  ED.  Fig.  4(c)  shows  the 
stress-strain  response  when  the  plate  is  compressed  along  the  ND. 
The  main  differences  between  the  response  along  the  ND  and  the 
other  two  directions  are  the  lack  of  a  well-defined  yield,  a  con- 
cave-down,  instead  of  an  s-shaped  hardening  curve,  and  greater 
strain  rate  sensitivity  in  both  the  initial  yield  strength  and  sub¬ 
sequent  hardening  behavior.  But  the  work  hardening  rates  are 
similar  to  the  TD  behavior.  The  yield  strengths  for  quasi-static  and 
dynamic  strain  rates  are  approximately  250  and  280  MPa, 
respectively. 

The  quasi-static  and  dynamic  mechanical  behavior  of  the  bar 
geometry  is  shown  in  Fig.  5.  As  in  the  previous  figure,  the  model 
simulations  are  plotted  but  will  be  discussed  in  the  next  section. 
The  mechanical  behavior  in  the  ED  has  the  same  nominal  s-shaped 
curve  as  the  plate,  with  a  yield  strength  of  approximately  340  MPa 
and  a  perfectly  plastic  stress  plateau  until  about  5.5%  true  strain. 
The  difference  in  the  work  hardening  rate  between  the  strain  rates 
is  more  pronounced  in  the  bar  as  compared  to  the  plate.  Ad¬ 
ditionally,  the  strain-to-failure  is  higher  in  the  quasi-static  strain 
rate  versus  the  dynamic  rate.  Loading  in  the  TD  (Fig.  5(b))  results 
in  a  plateau  that  is  barely  perceptible,  but  the  yield  strength  is  still 
the  same  between  the  strain  rates  (300  MPa).  The  work  hardening 
is  less  for  the  TD  of  the  bar  than  for  the  plate.  Finally,  the  behavior 
of  the  bar  in  the  ND,  Fig.  5(c),  shows  a  more  pronounced  plateau 
region  than  for  the  bar  TD  loading  direction;  this  curve  most  re¬ 
sembles  the  TD  of  the  plate.  Both  strain  rates  have  the  same  yield 
point  (290  MPa)  and  show  a  perfectly  plastic  flow  response  until 
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Fig.  3.  Pole  figures  with  randomly  sampled  orientations  of  AMX602  from  the  (a)  plate  and  (b)  bar  used  as  input  for  the  reduced-order  HCP  crystal  plasticity  model. 


approximately  3.5%  true  strain.  The  ensuing  work  hardening  and 
the  dependence  on  the  strain  rate  are  very  similar  to  the  TD  of  the 
bar.  The  strain-to-failure  continues  the  same  trend  with  a  higher 
value  at  the  quasi-static  strain  rate. 

We  want  to  highlight  three  overall  trends.  First,  the  difference 
in  the  yield  strengths  along  the  three  directions  for  each  geometry 
are  relatively  low  compared  with  conventionally  processed,  as 
well  as  ECAPed  Mg  alloys,  despite  sample  directions  coinciding 
with  peak  texture  components.  Two  possible  explanations  for  the 
reduced  anisotropy  are  that  a  weak  texture  is  formed  during 
processing,  and  that  the  processing  produces  a  relatively  small 
grain  size.  Although  the  preferential  textures  can  be  used  to  ra¬ 
tionalize  the  difference  in  responses  along  the  three  directions,  the 
reduction  in  intensity  compared  with  conventionally  processed 
Mg  alloys  implies  the  difference  in  responses  should  be  reduced. 
Additionally,  several  authors  [8,9]  have  postulated  that  Hall-Petch 
strengthening  is  greater  for  twinning  than  for  slip-based  me¬ 
chanisms.  Therefore,  as  grain  size  is  reduced,  so  is  the  relative 
difference  in  strengths  between  extension  twinning  and  slip- 
based  mechanisms. 

The  second  trend  is  between  the  geometries:  the  yield  strength 
and  ductility  at  the  quasi-static  strain  rates  are  greater  in  the  bar 
than  the  plate  samples.  The  former  is  because  the  extrusion 
temperature  of  the  bar  was  lower  (250°  versus  370°  C),  which 


locks  in  the  material  at  a  higher  energy  state  or  dislocation  content 
because  it  allows  less  recovery  of  the  microstructure.  The  latter 
may  be  caused  by  greater  porosity  or  powder  that  was  not  fully 
sintered  when  the  plate  was  manufactured,  as  the  larger  dimen¬ 
sions  permit  larger  temperature  and  stress  gradients  to  arise. 
Additionally,  the  extrusion  pressure  and  extrusion  ratio  (ratio  of 
initial  billet  cross-sectional  area  to  final  cross-sectional  area  fol¬ 
lowing  extrusion)  were  important  parameters  that  were  different. 

The  final  trend  is  there  is  good  agreement  between  the  bar 
geometry  here  and  the  results  of  Shen  et  al.  [30]  and  Elsayed  et  al. 
[27]  (the  bars  are  from  the  same  batch).  The  quasi-static  behavior 
is  similar:  their  results  labeled  “transverse  direction”  more  closely 
matches  our  ND  behavior.  But  they  meant  “transverse”  as  in 
transverse  to  the  ED,  which  is  true  for  both  the  TD  and  ND.  Their 
yield  strengths  at  the  dynamic  loading  rate  were  greater  than  the 
quasi-static  rate,  while  our  data  shows  no  difference.  An  obvious 
reason  is  their  strain  rate  was  4000 s_1  versus  1000  s_1  in  this 
study.  We  discuss  the  implications  of  the  different  strain  rates  in 
more  detail  in  Section  4.4. 

4.3.  Model  predictions  of  AMX602 

Simulations  were  carried  out  on  a  cylindrical  specimen  with  a 
height  to  diameter  ratio  of  1 :2  to  prevent  any  buckling  effects.  The 
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Fig.  4.  Mechanical  behavior  and  modeling  of  AMX602  plate  in  the  (a)  ED,  (b)  TD  and  (c)  ND. 


cylinder  was  discretized  using  approximately  1400  elements  and 
contained  approximately  300  grains.  Increased  mesh  refinement 
and  number  of  grains  was  found  to  have  little  effect  of  the  simu¬ 
lated  macroscopic  mechanical  response.  Representative  poly¬ 
crystalline  microstructures  were  instantiated  via  Voronoi  tessel¬ 
lation  with  grain  orientations  randomly  sampled  from  the  ex¬ 
perimentally  measured  pole  figures.  The  nominal  strain  rate  dur¬ 
ing  simulations  was  identical  to  those  provided  experimentally. 
Strain  to  failure  is  not  captured  in  the  simulations  as  the  post¬ 
bifurcation,  strain  localization  behavior  is  not  explicitly  modeled  in 
this  work.  Figs.  4  and  5  give  the  simulated  macroscopic  mechan¬ 
ical  response  of  the  polycrystalline  aggregates  that  represent  the 
AMX602  plate  and  bar  geometries  subjected  to  quasi-static  and 
dynamic  compressive  strain  rates.  There  are  several  noticeable 
similarities  and  differences  between  the  simulated  and  measured 
responses  that  require  further  discussion. 

The  most  pronounced  difference  is  that  the  model  is  unable  to 
capture  the  abrupt  yield  point,  particularly  when  the  material  is 
thought  to  yield  via  extension  twinning.  The  model  is  unable  to 
capture  such  an  abrupt  yield  point  because  extension  twinning 
and  basal  slip  alone  cannot  satisfy  the  von  Mises  compatibility 
criterion  required  to  accommodate  arbitrary  deformation  and 
cause  macroscopic  yielding;  however,  clearly  the  material  is  able 
to  yield.  Two  possibilities  for  these  differences  are  inherent  in 
model  approximations,  most  of  which  are  common  to  other  crystal 
models  of  HCP  metals  found  in  literature.  First,  the  model  assumes 
slip  confined  to  their  respective  glide  planes.  If  dislocations  were 
allowed  to  cross  slip  or  form  glissile  partials  that  act  as  dislocation 
sources  on  intersecting  slip  systems,  glide  on  a  single  system 
would  produce  multiple  independent  deformation  modes  and 
cause  yielding  with  fewer  slip  systems.  Second,  the  model  assumes 
that  deformation  due  to  twinning  can  be  approximated  in  a 


“pseudo-slip”  framework  that  smears  out  coordination  between 
finite  twinned  domains.  In  reality  twins  appear  to  form  co¬ 
ordinated  patterns  that  demarcate  finite  domains  of  twinned  and 
untwinned  material.  Direct  treatment  of  coordinated  deformation 
via  sequential  lamination  has  been  shown  to  provide  additional 
deformation  modes  and  energy  relaxation  as  compared  to  the 
averaged  response  [38,39  .  As  a  practical  substitute  to  soften  the 
material  response,  and  therefore  capture  the  abrupt  yield  behavior 
and  subsequent  plateau  when  the  material  deforms  via  extension 
twinning,  some  authors  40,41]  have  employed  various  relaxed- 
constraints  methods  [42,43]  with  relative  success.  Although  such 
methods  are  able  to  reproduce  initial  yield  behavior  while  em¬ 
ploying  the  aforementioned  restrictive  kinematic  approximations, 
they  act  as  a  substitute  to  the  rich,  complex  kinematics  exhibited 
by  real  materials,  and  do  not  solve  our  deficiency  in  addressing 
them  via  existing  continuum  crystal  descriptions.  Therefore,  we 
choose  to  use  the  more  accurate  kinematic  description,  i.e„  direct 
numerical  simulation  of  polycrystalline  aggregates. 

Another  difference  is  that  for  yield  behavior,  the  model  does  not 
predict  a  significant  difference  in  the  yield  magnitudes  between 
the  plate  and  bar  geometries,  whereas  the  bar  generally  has  higher 
initial  strength  in  each  direction.  The  model  parameters  were  fit  to 
the  plate  response,  and  no  difference  in  processing  was  included 
when  identifying  material  parameters.  The  model  is  only  simulat¬ 
ing  the  difference  in  the  textures  between  the  two  geometries,  but 
due  to  differences  in  the  processing  parameters  between  them, 
non-textural  effects  are  manifested  in  the  experiments.  One  dif¬ 
ference  is  the  extrusion  temperature  for  the  bar  was  120  °C  lower 
than  the  plate,  resulting  in  a  grain  size  which  is  1-2  pm  smaller, 
which  would  increase  the  yield  strength.  The  temperature  may  also 
modify  the  properties  of  the  precipitates  and  dislocation  structures 
within  the  material,  which  this  model  does  not  capture. 
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True  Strain 

Fig.  5.  Mechanical  behavior  and  modeling  of  AMX602  bar  in  the  (a)  ED,  (b)  TD  and  (c)  ND. 


Unlike  initial  yield,  the  work  hardening  behavior  is  captured 
well  by  the  model.  The  work  hardening  rate  in  the  ED  increases 
and  is  high  following  the  plateau  in  the  flow  stress  between  strains 
of  approximately  5-10%.  It  then  decreases  (but  is  still  positive) 
above  about  11%  until  failure  of  the  specimen.  Absolute  differences 
of  the  hardening  between  the  experiments  and  simulations  could 
be  due  to  the  inability  to  reproduce  the  initial  plateau  region.  The 
model  is  very  successful  at  capturing  the  experimental  hardening 
in  the  TD  of  the  bar  and  plate  and  in  the  ND  of  the  bar.  Even  the 
differences  in  the  stress  values  are  small.  The  deviations  are  larger 
when  loading  in  the  ND  of  the  plate  where  the  curves  are  more 
typical  concave-down  in  shape,  i.e.  slip  dominated  mechanisms. 
Finally,  the  model  captures  well  the  essential  differences  in  the 
mechanical  behavior  due  to  the  texture  difference  between  the 
plate  and  bar.  In  the  TD,  the  model  captures  well  the  minimal 
plateau  in  the  bar  and  a  more  pronounced  one  in  the  plate.  This  is 
similarly  true  for  the  ND  but  the  plate  lacks  a  plateau  region  while 
the  bar  has  one. 

With  reasonable  agreement  between  the  predicted  and  mea¬ 
sured  mechanical  response  of  AMX602,  particularly  in  the  post¬ 
yield  regime,  the  model  can  provide  further  insight  into  the  un¬ 
derlying  mechanisms  during  plastic  deformation.  Fig.  6  plots  the 
relative  activity  of  each  mechanism  for  the  different  loading  con¬ 
ditions  in  the  plate  and  bar.  Fig.  6(a),  (b)  and  (c)  corresponds  to 
loading  in  the  ED,  TD  and  ND,  respectively.  The  abrupt  onset  of 
twinning  causes  the  sudden  yielding  observed  in  Figs.  4(a)  and  5 
(a);  however,  the  twinning  contribution  quickly  reaches  a  max¬ 
imum  after  a  few  percent  true  strain.  Here,  the  mechanical  re¬ 
sponse  exhibits  an  inflection  point  as  twinning  exhausts  and  can 
no  longer  accommodate  the  applied  deformation,  giving  way  to 


slip.  Additionally,  the  model  and  texture  are  able  to  explain  why 
the  inflection  point  in  the  work  hardening  in  the  TD  occurs  at  a 
lower  plastic  strain  than  in  the  ED.  The  model  indicates  that  the 
contribution  of  twinning  is  greater  in  the  ED,  and  the  texture  is  the 
underlying  reason.  Twinning  is  greater  when  loading  in  the  ED 
because  a  greater  number  of  basal  planes  are  perpendicular  to  the 
ED  than  the  TD.  Twinning  accommodates  the  c-axis  extension 
when  the  material  is  compressed  perpendicular  to  the  c-axis. 
Extension  twinning  contributes  relatively  little  to  the  overall  de¬ 
formation  when  the  material  is  compressed  along  the  ND  (Fig.  6 
(c)).  Extension  twin  activity  is  suppressed,  since  for  the  idealized 
texture,  the  c-axis  undergoes  compression.  In  this  case  deforma¬ 
tion  is  primarily  accommodated  by  slip  mechanisms,  which  pro¬ 
duce  a  conventional  parabolic,  instead  of  s-shaped  hardening 
curve. 

In  each  loading  direction  the  plate  (rolling  texture)  versus  bar 
(extrusion  texture)  geometries  results  in  distinct  effects  on  the 
deformation  mechanics.  The  model  predicts  that  non-basal  slip  is 
the  dominant  mechanism  in  the  ED  at  the  onset  of  plasticity  (i.e.  at 
yield),  however,  due  to  its  high  initial  hardening  rate  it  gives  way 
to  basal  slip  and  twinning.  The  texture  difference  is  predicted  to 
have  only  a  small  effect  when  the  material  is  compressed  along 
the  ED  (Fig.  6(a))  because  for  both  cases  the  applied  load  is  per¬ 
pendicular  to  the  c-axis,  which  is  consistent  with  the  data.  The 
strain  rate  is  predicted  to  have  a  greater  effect,  but  this  is  not 
observed  in  the  data.  Where  the  non-basal  slip  is  greater  (bar),  the 
model  predicts  greater  flow  stress  differences  due  to  the  strain 
rate  sensitivity  of  these  mechanisms,  which  matches  the  data. 
When  specimens  are  compressed  along  the  TD,  the  model  predicts 
large  differences  in  the  amount  of  twinning  between  the  plate  and 
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True  Strain 


True  Strain 


True  Strain 


Fig.  6.  Relative  activity  of  deformation  mechanisms  in  the  (a)  ED,  (b)  TD  and  (c)  ND  as  predicted  by  the  model. 


True  Strain 

Fig.  7.  Mechanical  behavior  of  the  bar  in  the  ED  showing  the  strain  rate  in¬ 
sensitivity  of  extension  twinning. 

bar  (Fig.  6(b)).  In  the  bar,  more  c-axes  are  orientated  parallel  to  the 
loading  direction  than  in  the  plate.  Therefore,  slip  mechanisms  are 
enhanced  in  the  bar,  causing  the  s-shape  to  become  almost  non¬ 
existent.  Similarly,  the  rate-dependent  non-basal  mechanisms  are 
more  active  in  the  bar  than  plate,  resulting  in  the  greater  strain 
rate  sensitivity  observed  in  the  bar.  When  plate  specimens  are 
compressed  along  the  ND  (Fig.  6(c)),  the  material  exhibits  the 
highest  relative  activity  of  slip,  and  has  a  typical  concave  down 
stress-strain  response.  Specimens  taken  from  the  bar  exhibit  re¬ 
latively  more  twinning  and  an  associated  s-shaped  hardening 
curve,  similar  to  loading  in  the  plate  TD,  where  the  amount  of 
twinning  is  essentially  the  same  and  gives  nearly  an  identical 
mechanical  response.  The  model  predicts  greater  strain  rate  sen¬ 
sitivity  in  the  ND  direction  of  the  plate  because  the  non-basal  slip 
has  the  highest  relative  activity  at  true  strain  above  several 


percent  versus  the  other  two  directions,  which  is  reflected  in  the 
data.  Conversely  for  the  bar,  the  ED  data  show  the  greatest  rate 
sensitivity  but  the  TD  is  predicted  to  have  the  highest  relative 
activity  of  rate-sensitive,  non-basal  slip  above  several  percent 
strain.  Overall,  the  model  predicts  relative  activity  (or  the  trends  in 
the  activity)  of  the  different  deformation  mechanisms  between  the 
geometries,  loading  directions,  and  strain  rates,  in  a  manner  that  is 
consistent  with  the  measured  mechanical  responses  in  this  study. 

4.4.  Twinning  strain  rate  sensitivity 

In  all  cases,  except  for  the  plate  ND  behavior,  the  quasi-static 
and  dynamic  yield  points  were  essentially  identical  despite  several 
orders  of  magnitude  difference  in  the  strain  rate.  This  shows  that 
twinning  is  rate  insensitive  over  the  tested  strain  rates,  which  the 
model  assumes  (the  plate  ND  yield  point  is  slip  dominated).  Shen 
et  al.  [30]  shows  higher  yield  strengths  at  dynamic  versus  quasi¬ 
static  strain  rates  in  all  three  loading  directions  for  the  bar.  The 
difference  is  our  dynamic  experiments  were  at  1000  s_1  while 
theirs  were  4100  s-1.  We  performed  higher  rate  experiments  to 
see  if  extension  twinning  becomes  strain  rate  sensitive  at  a 
threshold  strain  rate  between  1000  s~]  and  4100  s_1.  Fig.  7  shows 
stress-strain  curves  at  different  strain  rates  in  the  ED  of  the  bar 
and  plots  the  corresponding  data  from  Shen  et  al..  At  3200  s_1  (the 
curve  is  not  shown  for  clarity)  there  was  no  indication  of  twinning 
becoming  rate  sensitive,  where  it  very  closely  matched  the 
1000  s~]  curve.  At  5000  s~]  there  is  evidence  for  an  increase  in 
the  abrupt  yield  point  (a  3800  s_1  experiment  is  very  similar  to 
the  5000  s_1  but  is  not  shown  for  clarity).  The  yield  point  matches 
that  from  Shen  et  al.  but  the  stress  drops  down  to  a  magnitude  the 
same  as  at  lower  strain  rates  (at  approximately  3.5%  true  strain), 
whereas  the  post-yield  behavior  is  perfectly  plastic  for  Shen  et  al.. 
Initially  we  concluded  this  was  evidence  for  extension  twinning  in 
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this  material  becoming  strain  rate  sensitive  between  3200  s  1  and 
3800  s-1. 

But  the  observed  “hump”  around  yielding  and  the  subtle  os¬ 
cillation  in  the  vicinity  of  the  arrow  could  be  wave  dispersion  ef¬ 
fects  between  the  high  impedance  bars  and  the  lower  impedance 
Mg  sample,  which  can  cause  oscillations  in  the  stress-strain  curve 
[32],  Or  error  may  be  generated  because  the  magnitude  of  the 
transmitted  wave  (which  is  used  to  calculate  stress)  was  small 
relative  the  magnitude  of  the  reflected  wave  (which  is  used  to 
calculate  strain  and  strain  rate).  The  “hump”  and  subtle  oscillations 
were  present  at  3800  s  1  also.  We  attempted  to  reduce  this  error 
by  testing  a  sample  with  a  reduced  (loading  direction)  thickness. 
The  strain  and  strain  rate  are  inversely  proportional  to  the  sample 
thickness,  so  to  achieve  the  same  strain  rate  as  a  thicker  sample, 
less  input  stress  needs  to  be  generated  in  the  incident  bar,  i.e.  the 
incident  wave  magnitude  is  reduced.  The  transmitted  wave  mag¬ 
nitude  stays  the  same  because  the  cross-sectional  area  was  not 
changed  but  the  reflected  wave  reduces  in  magnitude,  which  im¬ 
proves  their  relative  magnitudes.  A  representative  experiment 
with  a  reduced-thickness  sample  is  shown  at  a  strain  rate  of 
6000  s_1  and  the  yield  is  essentially  the  same  as  at  1000  s  \  i.e. 
no  “hump”,  and  there  is  no  subtle  oscillations  in  the  subsequent 
work  hardening  regime  of  the  curve.  A  few  experiments  with  re¬ 
duced-thickness  samples  at  strain  rates  below  6000  s  1  are  con¬ 
sistent  with  the  curve  shown.  We  should  note  that  our  data  has 
not  been  filtered  or  smoothed  in  any  way-the  signals  were  only 
multiplied  by  the  requisite  constants  to  obtain  stress,  strain  rate 
and  strain. 

Thus,  our  final  conclusion  is  extension  twinning  for  this  Mg 
alloy  is  most  likely  not  strain  rate  sensitive  between  the  strains 
rates  tested  here  (10-4  to  ~5000s-1),  which  supports  the  as¬ 
sumption  of  this  model.  The  data  from  Shen  et  al.  hints  there  could 
be  some  sensitivity,  but  the  yield  strength  increase  is  less  than 
10%,  so  it  could  be  variability  in  the  behavior  between  samples  or 
wave  dispersion  at  low  plastic  strain  producing  anomalous  results. 
It  is  possible  twinning  does  become  rate  sensitive  to  a  greater 
effect  at  rates  beyond  those  tested  here,  but  the  effects  of  the 
material,  grain  size,  texture,  anisotropy,  etc.  on  the  degree  of 
sensitivity  are  unknown  and  require  further  study. 


5.  Conclusions 

In  this  study  a  high  strength  Mg  alloy,  AMX602  (Mg-6%Al-0.5% 
Mn-2%Ca),  was  produced  by  the  spinning  water  atomization  pro¬ 
cess  (SWAP)  and  subsequently  extruded  into  a  plate  or  bar.  The 
microstructure  was  fine-grained  (the  vast  majority  of  grains  be¬ 
tween  0.5  and  5  pm)  with  relatively  weak  texture  versus  other  Mg 
alloys.  The  texture  differed,  from  an  extrusion  texture  in  the  bar  to 
a  rolling  texture  in  the  plate.  Quasi-static  and  dynamic  compres¬ 
sive  loading  showed  AMX602  possesses  high  strength,  reasonable 
ductility  and  low  anisotropy  in  the  principal  processing  directions. 
Crystal  plasticity  modeling  of  the  rate-dependent  mechanical  re¬ 
sponse  was  in  relative  agreement  with  the  measured  experimental 
data,  especially  when  all  the  deformation  mechanisms  are  active 
in  abundance  simultaneously.  The  model  was  used  to  parse  the 
relative  contribution  of  each  deformation  mechanism  for  given 
loaded  cases,  and  study  how  texture  differences  between  different 
extruded  geometries  would  affect  quasi-static  and  dynamic  mac¬ 
roscale  behavior.  In  the  plate,  the  model  indicated  extension 
twinning  was  profuse  for  compression  along  the  ED  and  TD, 
whereas  slip  dominated  the  response  for  compression  along  the 
ND.  The  mechanisms  of  the  bar  are  similar  to  the  plate,  but  in  the 
ND  it  twins  enough  to  have  the  same  behavior  as  the  TD.  Finally, 
we  see  little  to  no  evidence  for  extension  twinning  becoming 
strain  rate  sensitive  over  the  range  of  strain  rates  tested  here, 


which  supports  our  model  assumption. 
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